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The spreading of an epidemic is determined by the connectiviy patterns which underhe the pop- 
ulation. While it has been noted that a virus spreads more easily on a network in which global 
distances are small, it remains a great challenge to find approaches that unravel the precise role 
of local interconnectedness. Such topological properties enter very naturally in the framework of 
our two-timestep description, also providing a novel approach to tract a probabilistic system. The 
method is elaborated for SIS-type epidemic processes, leading to a quantitative interpretation of the 
role of loops up to length 4 in the onset of an epidemic. 
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I. INTRODUCTION 

Almost all of us have already met someone far from 
home, who turns out to share a friend with us. Milgram 
first put this social phenomenon on a firm basis, finding 
that any two individuals are separated, on average, by six 
acquaintances 0. Furthermore, our society is composed 
of groups within which, if individual A is acquainted with 
persons B and C, then B is likely to know C. In the lan- 
guage of social networks, this means that the individu- 
als A, B and C are arranged on a triangle. These two 
facts tell us that in the social universe, global distances 
are small, and local interconnectedness, i.e. clustering, is 
high. A network that possesses both of these topologi- 
cal properties is called a small world. In order to find 
a model that accounts for both of these properties, one 
could choose a regular lattice whose nodes are indeed lo- 
cally highly interconnected. But since global distances 
are large in this type of network, it is an ineligible can- 
didate for a small world model. In a random graph ^j, 
that is a set of nodes (e.g. laid out on a virtual circle) 
with connections between them established at random, 
a link is more likely to point to a far away node than 
to one close by. (Here we refer to the distance on the 
virtual circle's circumference.) As a consequence, the 
presence of the very many long-range links account for 
the desired global property, but the degree of clustering 
is low. Watts and Strogatz combined these two insights 
and proposed a model that interpolates between a regu- 
lar lattice and a random graph, thus capturing both the 
global and local topological properties mentioned above 
Q ■ But the small world property is not merely exhibited 
by social networks. Yet, through the increased availabil- 
ity of data, it was found that the simultaneous occurence 
of high local and global interconnectedness is prominent 
to a much wider class of systems, notably the Internet 
H01, the World Wide Web metabolic networks jj 
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and food webs 

In addition to the small world property, there is an- 
other important fact when it comes to characterizing the 
topology of a complex network: not all the nodes have 
the same number of edges. The corresponding measure 
is the degree distribution P{k) which gives the probabil- 
ity that a randomly chosen node has degree k, that is 
k edges. For most of the just mentioned examples, this 
distribution was found to be P{k) ^ k^^ , 2 < 7 < 3, 
implying the absence of a characteristic (degree) scale, 
hence the name scale-free network. This emergence of 
scaling can be understood, for example, in terms of grow- 
ing networks. Starting from a small core graph, at each 
time step a node is added together with a certain number 
of edges that are connected to existing nodes, the latter 
being chosen preferentially, that is a link is more likely 
established to a high degree node |3- This fat-tailed de- 
gree distribution implies the presence of hubs (high degree 
nodes) which hold together the network and play a cru- 
cial role in issues such as robustness or fragility [lOl . 

This explosion of research activity in the field of com- 
plex networks has also shed new light on epidemic spread- 
ing since the latter can be regarded as a dynamical pro- 
cess occuring on a complex network: a computer virus 
spreads on the Internet, and HIV (as a biological ex- 
ample) propagates on top of the web of human sexual 
contacts. Also these two examples fall into the category 
of scale- free networks 0|, and it is this topologi- 
cal property that accounts for the absence of a finite 
epidemic threshold in the corresponding spreading phe- 
nomenon p^ . 

The degree distribution is only a first way to charac- 
terize the degree related topology of a complex network. 
Indeed, by analyzing scientific collaboration networks, re- 
searchers working with many others (high degree nodes) 
tend to collaborate with other "hubs" . This means that 
there exist degree correlations, and the just mentioned 
property has been called assortative mixing, holding gen- 
erally for social networks ,15.] . On the other hand, in the 
Internet (at the autonomous system level), high degree 
nodes are more likely connected to low degree ones, thus 
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exhibiting disassortative mixing. The influence of such 
degree correlations on the spreading of an epidemic was 
investigated in detail finding that neither assorta- 
tive nor disassortative scale-free networks exhibit a finite 
epidemic threshold 0|. The insights about the role of 
these degree related topological properties in epidemic 
spreading have been gained at the mean-field level. 

Besides degree correlations, triangles are ubiquitous in 
complex networks as outlined in the first paragraph, and 
more generally, many loops of short length were found 
in these systems Motivated by Watts and Strogatz' 
model 0] which uses a regular lattice, i.e. an ordered net- 
work possessing many loops, as starting point, we shall 
also use the concept of local ordering when referring to 
the loop structure of a complex network. More informa- 
tion can be extracted if the triangles are sorted according 
to the degrees of their corners, finding that mainly the 
low deg ree nodes account for the high level of clustering 
[l^ l2(ll |. This suggests the presence of interesting modu- 
lar organizations, and similar results were obtained from 
an analysis of loops of length 4 ^| . Local ordering plays 
a crucial role in the function of a metabolic network (with 
scale- free topology). Indeed, the loop structure is much 
richer than in the scale-free model based on growth and 
preferential attachment This work also investigates 
the number of triangles as a function of the system size 
for this model, a result which was generalized to loops 
up to length 5, yielding robust scaling relations [2^ . 

Obviously the presence of loops has an effect on the 
spreading behavior since, with respect to a treelike topol- 
ogy, there exist many more paths along which the virus 
can propagate. Different strategies in order to gain in- 
sights about the role of local ordering properties have 
been proposed. An interpretation of how clustering in- 
fluences the stationary spreading behavior was obtained 
by mapping the epidemic process onto bond percolation 
[23 . Another approach is to abandon the mean- field level 
and take into account spatial correlations which govern 
the epidemic dynamics. Matsuda et al. first used the or- 
dinary pair approximation in order to study a population 
dynamical problem . This approximation, as its name 
anticipates, accounts for pair correlations and lies at the 
basis of improved pair models [25I |2^ which uncover the 
role of local ordering in a rather indirect way: cluster- 
ing enters by making a number of assumptions about the 
open (Z) and closed (A) triple correlations. In cluster 
approximations, a time-dependent probability is assigned 
to each configuration of the fundamental cluster whose 
choice is guided by the network topology [23| : for investi- 
gating the spreading dynamics on a triangular lattice one 
uses a triangle as fundamental cluster whereas the star 
is the appropriate choice for a random network. Higher- 
order correlations are therefore embedded very explicitly. 
Moreover the systematic improvability of this method 
makes it a powerful tool to study probabilistic systems. 

In order to understand the role of local ordering prop- 
erties, the exploration of temporal correlations seems to 
be an even more natural approach: for example within a 



two-step description, it matters whether the local topol- 
ogy is treelike or if loops of short length are present. 
The method is illustrated for the susceptible-infected- 
susceptible model (see for example 28]), homogeneous 
networks are used as starting point, and analytical es- 
timates are obtained also for disordered graphs obeying 
P{k) = Sk,K, K being arbitrary. 

The paper is organized as follows. Section II describes 
the adopted model consisting of the contact network as 
well as the local dynamics. In section III, we introduce 
the formalism, from which the two-step description is 
derived. For completeness, the necessary ingredients in 
order to arrive at the one-step site approximation, i.e. 
the common mean-field level, are also shown. Section 
IV explores the implications of the double-step approach 
for networks of degree 4, the generalization to arbitrary 
degree is done in the following section. The major con- 
clusions are drawn in section VI. 



II. THE MODEL INGREDIENTS 

The dynamical laws that describe the spreading of an 
infectious disease are determined by the contact struc- 
ture which underlies the population. We therefore model 
the epidemic as a dynamical process on top of a given 
network that does not change in time. The nodes of the 
network represent individuals, and the links correspond 
to relationships between individuals along which an in- 
fective agent can propagate. 

Since the aim of this paper is the investigation of the 
role of loops of short length, we adopt a rather simple 
epidemiological model where the individuals can be only 
in two possible states, namely infected (I) or susceptible 
(S). Because the nodes repeatedly run through the cy- 
cle susceptible — > infected susceptible, it is called SIS 
model. In the physics community, it has recently been 
formulated as follows Ti^ : A node susceptible to the dis- 
ease gets infected with probability i^At if it is connected 
to at least one infected nearest neighbor. On the other 
hand, infected nodes recover spontaneously with prob- 
ability 5/S.t. This version of the SIS model is formally 
advantageous with respect to its conventional formula- 
tion, where infected nodes can infect neighboring sus- 
ceptible vertices with probability v/Si In the latter 
case, susceptible nodes become infected with probabil- 
ity 1 — (1 — ;yAt)'^'°f, fcinf being the number of infected 
nearest neighbors. In this paper, we will use the former 
version. By rescaling the time unit, we can reduce the 
number of parameters to one: the time evolution is de- 
termined by the effective spreading rate A = v/S, and 
the recovery rate is set to 1. The quantitative details of 
the behavior of the system still depend on the choice of 
At. In particular, the effect of the loops are of higher 
order in At, such that their influence is not seen in the 
continuous-time limit (At ^0). As long as At > 0, we 
set this quantity to 1 without lack of generality. 

The other model constituent concerns the underlying 
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(a) (b) (c) (d) 

FIG. 1: Homogeneous networks of degree 4. In (a) the nodes are connected at random under the restriction that K — 4 links 
emanate from every vertex, leading to a treelike topology. The role of triangles and loops of length four is studied by means of 
the square lattice (b) and the Kagome lattice (c) where they appear separately, as well as the ring-type network (d). 



network. We shall not attempt to examine the combined 
effect of the degree distribution, degree correlations and 
the loop structure. While the role of these degree related 
connectivity patterns in epidemic spreading is rather well 
understood, little attention has been given to the effect 
of local ordering. That is why we focus on this topolog- 
ical property. We therefore use networks in which every 
node has the same number of nearest neighbors (fixed 
degree). The adopted strategy is to start with strictly 
homogeneous graphs (networks in which identical con- 
nectivity patterns are "seen" from every node) differing 
in the detailed loop structure (Fig.^. This leads us to an 
understanding of the role of these distinct local ordering 
properties, even for networks that are nolonger strictly 
homogeneous but which still obey P{k) — 6k, k, K being 
the constant degree. 

In the case K = A, examples of purely homogeneous 
networks are the random homogeneous network (Fig. 
ipi) constructed according to the MoUoy Reed algorithm 
|29', '3^, the square lattice (Fig.^), the Kagome lattice 
(usually used in condensed matter physics as this geom- 
etry represents one of the most frustrated antiferromag- 
netic systems, Fig. ^) and a ring where nodes two units 
apart from each other are also directly connected (Fig. 
^). For the random network, the homogeneity defini- 
tion given above only holds approximately. While loops 
of short length do not occur in the limit of large network 
size, long ones of various lengths may exist. Therefore 
different nodes may not "see" identical connectiviy pat- 
terns. Also the Cayley tree lacks completely in loops. 
But as the last generation of nodes has degree 1, it does 
not fall into the class of homogeneous networks. Indeed, 
since the number of nodes belonging to generation N 
is K{K — most of the nodes are even comprised 

within the last generation. Clearly, all the networks de- 
picted in Fig. n are characterized by the degree distri- 
bution P{k) = 6k4, but they differ in the way the sec- 
ond neighbors are arranged, that is the loop structure 
becomes richer going from (a) to (d). While the random 
network has no short loops at all, the ones of length 4 are 



a fingerprint of the square lattice. In the Kagome lattice, 
every node is a corner of two distinct triangles. Finally, 
in the ring of degree 4, triangles and two different kinds of 
loops of length 4 are found, namely plaquettes (as in the 
square lattice) and the so-called primary quadrilaterals 
(two adjacent triangles), see section IVB. 

In order to study the epidemic process on disordered 
"homogeneous" topologies, we subject the networks in 
Fig. [TjD-d to the following rearrangement algorithm (Fig. 
m 31]: 

• Choose randomly two links (link 1 connecting node 
Ai with Bi and connection 2 linking vertex A2 with 
B2) that do not share a common node. 

• Remove these two links and establish two new con- 
nections between Ai and B2 as well as A2 and Bi. 

Repeating this rewiring procedure a certain number of 
times leads to locally varying numbers of loops while the 
degree distribution remains unaltered. 



A, ©- 




FIG. 2: Rewiring procedure not affecting the degree distri- 
bution: the end vertices of two arbitrarily chosen links are 
exchanged. 



III. THE FORMALISM 

In this section, we shall introduce the formalism that 
is used in order to describe the dynamics on top of the 
networks specified in the previous section. For complete- 
ness, we then derive the (single-step) mean-field approx- 
imation which could also be obtained heuristically. In 
the last paragraph, we derive a two-step description that 
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allows us to gain insight about the loop structure in epi- 
demic spreading. 

On an exact level, we shall describe the epidemic dy- 
namics by assigning a probabilty 'Pt(x) to each configu- 
ration X at every instant of time t. The vector x contains 
the states Xi of all the nodes i of the network, Xi being 
either (susceptible) or 1 (infected) . The system proba- 
bilities satisfy at any time t 

X 

and evolve in time according to 



(1) 



The transition matrix of the system Wy^x is obtained 
from the matrix Wy^^^^ which shall denote the probabil- 
ity that the state of the arbitrary site I changes from yi 
to Xi, through 



N 



1=1 



N being the total number of nodes in the network. The 
matrix elements representing the probabilites for the pos- 
sible events at the site / are given by our version of the 
SIS model, namely 



W, 



0^1 



= A 



1 n(i 



jnnl 



wl 



1 - A 



jimi 



or in a more compact form 



= 1 - XI 



\{2xi^l){l-yi)\l-\{{l 



jnnl 



The products in the above expressions have to be taken 
over all the nearest neighbors j of node I. The factor 
^~Yljnni(^~yj) ^ least one yj = 1 and otherwise. 

Before deriving our two-timestep description, we show 
how the conventional mean-field approximation is re- 
trieved through this formalism. At that level, the sites 
are considered independently from each other, and we 
write for the system probability 



N 



rti^) = llPt{xi), 



(2) 



i.e. the system is described by the single variable -Pt(l) 
[the probabihty of being susceptible is Pt{0) — I — Pt{^)]- 
Its dynamics is obtained from Eq. |^ by summing it over 
all possible configurations x, xq held fixed 

J2 Pt+i(x)^5]7',(y) J2 



where the node can be chosen in an arbitrary way. The 
left hand side of the above equation corresponds to the 
probability that node is in state xq at time t+1. With 
the ansatz the time evolution is 



Pt+i(l) = APt(0)[l-Pt(0) 

With Pt{l) = Pt and Pt{0) = 1 — Pt, by consequence, we 
obtain for small values of Pt 



Pt+i^\K{l~Pt)Pt. 



(3) 



The stationary-state condition (Pt+i = Pt) leads to a 
value for the epidemic threshold 



A. = 



K 



(4) 



Therefore all the networks depicted in Fig. ^ are treated 
identically at this level of description. The only topo- 
logical property that determines the critical value of the 
effective spreading rate A is the number of nearest neigh- 
bors K. 

Improvements upon the mean-field description can 
be obtained by taking into account spatial correlations. 
Based on the ordinary pair approximation , the pres- 
ence of triangles enters by establishing a number of hy- 
potheses about the open and closed triple correlations 
[2^|2g. This approach therefore embeds this topological 
property in a rather implicit way. In cluster approxima- 
tions [23, the local topology and its associated spatial 
correlations translate directly into the choice of the clus- 
ter and a set of probabilities for its possible configura- 
tions. 

Another strategy that serves to incorporate local or- 
dering properties is to take into account temporal corre- 
lations. Thus by performing two timesteps exactly, we 
expect that the way the second neighbors are arranged, 
enters very naturally into the description. For example, 
the cases where two nearest neighbors of an arbitrary 
node are also directly connected (presence of a triangle) , 
where they are linked via a second neighbor (giving rise 
to a loop of length 4) or where the only path goes through 
the original node (treelike structure), lead to different re- 
sults. In the remaining part of this section, we derive the 
general equation, special cases are then looked at within 
the following section. 

As outlined above, we iterate Eq. ^ once 

Pt+l(x) - ^[Wy^xE W._yP*_l(z) 



We now again pass to a site approximation. By summing 
the above equation over all possible configurations x, a;o 
held fixed, we get 



Pt+i(xo)=E[Pt-i(z) 

^ {yi}i=o,i, 



K 

-0 n 

(5) 
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FIG. 3: Equilibrium prevalences for the epidemic process on top of different networks characterized by P{k) = Ski- Left: 
Simulation results for the homogeneous random graph (★), the Kagome lattice (A), the square lattice (□) and the ring-type 
network (o). Right; The one-step site approximation (thin solid) ignores (local) ordering properties and yields Ac = 0.25 for 
any homogeneous network of degree 4. For the MoUoy Reed network (dotted), the Kagome lattice (solid), the square lattice 
(dashed) and the ring (dotted-dashed), different steady-state prevalences are obtained at the two-step level. 



again being an arbitrarily chosen node. Furthermore, 
the system probabihty Vt-i{z) is given by Eq. (|2Jl, and 
the nodes 1,2,..,K denote the nearest neighbors of the 
arbitrarily chosen node 0. As only z-states associated to 



the vertex 0, its nearest and second neighbors appear in 
the VF-factors, the sum over the ^-variables associated to 
nodes more than 2 hnks away from vertex 0, is carried 
out triviaUy. A tour de force calculation leads to 



K K 



Ql— 1 Qi — 1 0:2— Cti + l Qi — 1 

K K K K K 

+ ^'(E E E (/"i/-./c.3)t-i+ E E (/o/-J-. 

Ql — 1 Q2 — CKl + l Ct3 — + I ai— la2— Q^l+l 

K K K K K K 

^(^'^)'^(E E E (/ai/a2---/a/f )t-l + E E '■■ E (/o/qi /as ■ ■ -/aff- 1 ) t-1 



ai — 1 OL2—OLI + 1 aK—0£K-l+^ 

Thereby the connectivity embedding factor 

u={i~ [1 - n (1 " 

anna. 

J 1 if = and at least one Z(j — 1^ 
I otherwise. 



Ql — 1 a2— (ii + l a/^_i— aK-2+1 

I 



(6) 



node 0, its nearest and second neighbors, these vertices 
collectively being denoted by A/"^ , 

(7) 

were introduced for notational convenience. In the fol- 
and the expectation value of a function of the states of lowing, an expectation value of a product of n /-factors 
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FIG. 4: Pictorial representation of the terms of Eq. Q for a treelike topology. The order of a specific subgraph is given by the 
number of filled circles, (a) corresponding to (/i), (b) to (/1/2) and (/o/i), and so forth. In this vein, the indices of the /-factors 
appearing in the expectation values correspond to filled circles whereas empty circles represent their nearest neighbors. The 
dashed lines are the links of the complete graph which are not contained in a specific subgraph. 



will be referred to as a term of n-th order although it is 
proportional to A • A". As is illustrated in detail in the 
next section, every term of Eq. © corresponds to a sub- 
graph of the graph composed of the nodes N"^ and whose 
links are according to the network under investigation. It 
can already be anticipated that the first term accounts 
for the degree distribution only, whereas the contribu- 
tions of higher order will give insight about the role of 
the loop structure. 



IV. NETWORKS OF DEGREE 4 

In this section, we elaborate the implications of our 
two-step description, the analysis being restricted to the 
stationary state. 

In the left part of Fig. O we show the simulation re- 
sults for the graphs introduced in section II. The ring- 
type network exhibits the largest epidemic threshold. For 
both the square and Kagome lattices, the critical value 
is Ac ~ 0.34. We therefore anticipate that either four 
plaquettes or two triangles (per node) lead to the same 
effect in the regime of low prevalences. Finally, the lowest 
epidemic threshold is found if the population is arranged 
on a homogeneous random network (of degree 4). The 
last result is very intuitive since in such a graph, global 
distances are small, making it more easy for a virus to 
spread. Therefore, even if the effective spreading rate is 
rather low, a finite fraction of the population will be in- 
fected in the stationary state, hence the small value for 



the location of the onset of the epidemic. In summary, 
these results indicate that the poorer the loop structure, 
the lower the corresponding epidemic threshold. 

In the right part of Fig. O the one-step and two- 
step site approximations are reported. The former corre- 
sponds to the steady-state solution of Eq. ||2Jl for which 
p = 0atAc = l/4 according to Eq. igj. All the net- 
works in question are therefore treated identically, the 
loop structure being ignored at this level of description. 
Yet, the two-step solutions [Eq. |SJ|] are diverse for the 
different graphs. Going from right to left, the curves 
correspond to the ring, the square lattice, the Kagome 
lattice and the MoUoy-Reed network, that is they appear 
in the same sequence as at the level of simulation. Fur- 
thermore, the curves corresponding to the Kagome and 
square lattice also meet the oi-axis at the same value of A. 
These findings confirm our intuitive arguments given in 
the previous section. It has to be noted that the two-step 
estimates for the threshold values are still considerably 
inaccurate especially for the ring and lattices, but this 
just highlights the presence of higher-order spatiotempo- 
ral correlations. However, the important point is that the 
degeneracy associated to the one-step description disap- 
pers at the two-step level. 

On the basis of Eq. © , we shall now analytically study 
the effect of local ordering properties upon the epidemic 
spreading, leading to a quantitative understanding of the 
threshold value. 
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A. Random network 

We shall now evaluate all the terms of Eq. ® for a lo- 
cally treelike topology. Fig.^shows the subgraphs repre- 
senting the terms in Eq. ©, in increasing order. Thereby 
the correspondance is as follows: Given the term (fafp), 
the nodes a and /3 are represented by filled circles whereas 
their nearest neighbors are drawn by empty circles. The 
links which enter at the level of the subgraph in question, 
are represented by solid lines whereas the ignored ones 
are dashed. If we denote the second neighbors of the cen- 
tral vertex by 11,12,13 for I = 1,2,3,4 and follow Eq. 
(|7|). the first order contribution for ai = 1 (subgraph in 
Fig. Hi) is 

21 20 211. .213 

X (1 - Zl)[l - (1 - Zo)(l - Zll)(l - Z12)(1 - Z13)]| 

/l 

where the sum over the z-variables to which no circles are 
associated, has been carried out trivially. Furthermore, 
we again have set P = P{1) in the third line (this is 
also done below), and the time index was omitted since 
we are only interested in the steady state. This term 
appears with multiplicity 4 (due to X]^i=i)i gi'^iiig the 
contribution 16P to first order in P. 

Fig.^ shows the subgraphs representing (/1/2) (upper 
part) and (/o/i) (lower part). Their contributions are 

(/i/2)=EEE E E {^(^o)p(zi)p(^2) 

Zq Z\ Z2 Z11..Z13 Z21--Z23 

xP{zii)P{zi2)P{zi3)P{z2l)Piz22)P{z23)flf2} 

^P + OiP^), 
occuring (2) = 6 times and 

(/0/l)=EEE E {P(.Zo)Piz^)Piz2)P{z,)P{z,) 

Zq Z\ Z2--Z4^ ^Il-.^IS 

xP(zn)P(2i2)P(zi3)/o/i} 
= + C'(p3), 

thus not giving a contribution to first order in P. As 
a consequence, the second-order term (that is the one 
proportional to A • A^) is 6P. 

As the procedure should now be clear, we only give the 
results for the remaining orders. The upper subgraph of 
Fig. ^ represents the term {f 1/2 fa)- Its contribution is 
P-I-C'(P^). The lower subgraph corresponds to (/0/1/2), 
yielding 18P^ -I- C'(P^). As the former term has multi- 
plicity (3) ~ 4, the total third-order contribution is 4P. 

As far as the fourth order is concerned (Fig. 011) , the 
subgraph involving node as filled circle neither gives 





FIG. 5: Loops of length 4 in a network. Primary quadrilater- 
als (left) are two adjacent triangles and therefore involve only 
nearest neighbors (solid empty circles) of the central node 
(filled circle). A secondary quadrilateral (right) involves a 
second neighbor (dotted empty circle) as well. 



a contribution whereas the term (/1/2/3/4) having mul- 
tiplicity 1 also gives P -I- O(P^), thus totally yielding 
A • A^P. 

Collecting these findings, we obtain the following con- 
dition that determines the epidemic threshold for a tree- 
like topology 



1 = A(16A-6A2 + 4A3 - A^), 



(8) 



which is satisfied by Ac — 0.2609. This is the value that 
can be extracted from the right panel of Fig. O (second 
curve from the left). 



B. Graphs with loops 

It is easy to imagine that the preceeding analysis yields 
different results when triangles and loops of length 4 
are present. Caldarelli et al. have classified loops 
of length 4 in a complex network into primary and sec- 
ondary quadrilaterals (Fig. [SJ. In the former case, the 
external vertices which the loop is composed of are all 
nearest neighbors whereas secondary quadrilaterals are 
plaquettes, the external nodes being two nearest and a 
second neighbor. With these concepts, the loop struc- 
ture of a strictly homogeneous graph can quantitatively 
be characterized as follows: By choosing an arbitrary 
node, the number of edges between its nearest neighbors 
is denoted by E. Qi and Q2 shall refer to the number of 
primary and secondary quadrilaterals. For the networks 
depicted in Fig.^-d, we report the corresponding values 
in Tab. D 



TABLE I: Loop properties for the simple non-treelike net- 
works described in section II. 





E 


Qi 


Q2 


square lattice 








4 


Kagome lattice 


2 








ring 


3 


2 


2 



Let us now look at the subgraph development for the 
square lattice whereby we focus on the important changes 
with respect to the treelike case. The reader interested 
in the full elaborations is referred to the appendix. We 
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FIG. 6: Second-order subgraphs not involving the central 
node as filled circle for the square lattice. To first order in P, 
the left subgraph yields the contribution 2P, the right one P. 
For further explanations, see Fig.|l] 



have already noticed that the first-order term is fully de- 
termined by the degree distribution, therefore the A • A 
coefhcient is 16, as in the treelike case. At order 2, the 
term (/1/2) (upper subgraph of Fig. ^) splits into two 
subgraphs in the presence of plaquettes (Fig. 0). The 
right subgraph is the same as in the treelike case, yet the 
left yields a contribution 2P + 0{P^). Their multiphc- 
ities are 4 (left) and 2 (right) summing up to (2) — 6. 
The resulting A • A^ coefficient is therefore -10. Although 
different subgraphs enter into the development also at 
the orders > 3, the coefficients appearing in the equation 
determining the epidemic threshold do not change. 





FIG. 7: Subgraphs of second order for the Kagome lattice, 
both contribute with P + 0{P'^). See Fig. 21 for details re- 
garding line- and circle styles. 

The second-order subgraphs for the Kagome lattice are 
depicted in Fig.[7| Both's contributions are P + 0{P^). 
The one involving a triangle appears 6 times whereas the 
right subgraph has multiplicity 4. We therefore obtain 
the same second-order coefficient as for the square lat- 
tice. An analysis for the higher-order subgraphs yields 
no difference with respect to the square lattice. These 
two cases are therefore equivalent at the two-step level 
for P < 1. 

In Tab.im we summarize the coefficients for these two 
lattices as well as the ring (Fig. QJi) . The full develop- 
ments are given in the appendix. 

Of course the idea is now to extend Eq. © such that 
it holds for all the investigated graphs. Our findings sug- 
gest that the local ordering properties enter in the fol- 



TABLE II: Coefficients of the two-step threshold equation for 
our networks having in common P{k) = 5^4, but differing in 
the loop pattern. 



A ■ A"- coeff. 


n = 1 


n = 2 


n = 3 


n = 4 


square lattice 


16 


-10 


4 


-1 


Kagome lattice 


16 


-10 


4 


-1 


ring 


16 


-16 


6 


-1 



lowing way into the equation determining the epide 
threshold: 



;mic 



1 = A 16A-(6 + 2£; + gi + Q2)A2 + (4+Qi)A3-A'^ . (9) 






(a) 



(c) 



FIG. 8: Loops of length 5 involving different hierarchies of 
nearest neighbors. Connections emanating from the central 
node (filled circle) are drawn as a solid line, links going from 
nearest neighbors (empty solid circle) to other nearest neigh- 
bors or second neighbors (empty dotted circle) are dashed, 
and second neighbors are connected by a dotted line. The 
pentalateral (a) involves nearest neighbors only, in (b) the 
loop traverses a second neighbor and in the one in (c) lacking 
in internal connections, a link between two second neighbors 
serves to close it. 

However this is not the full story. What about loops 
of length 5? Let us argue why they do not enter in the 
framework of a two-step description. Although there ex- 
ist such loops involving only first and second neighbors 
(Figs.lHt and b), it may also be closed only between two 
second neighbors (Fig.lSt). Obviously such a connection 
is ignored at the two-step level. In the language of graph 
theory the latter case corresponds to a fundamen- 
tal loop whereas the former examples can be reduced to 
loops of length 3 and 4. But whatever the number of 
hierarchies of nearest neighbors involved in the forma- 
tion of the loop is, the point is the following. If the 
central node is infected at time t, it can causally affect 
only vertices two links away, corresponding to a chain of 
4 links. Obviously, it matters whether the first and the 
last node of this chain are identical. In this case, we have 
a loop of length 4. Otherwise it cannot be distinguished 
whether the topology is fully treelike of if loops of length 
greater than 4 are present. Along these lines, it has to be 
expected that loops up to length 2n enter within an n- 
timestep description. In contrast, the presence of higher- 
order quadrilaterals modifies the coefficients of Eq. 



fa) 



FIG. 9: Quadrilaterals of orders 3 and 4. The n (a: n = 3, b: 
n = 4) nearest neighbors (empty solid circles) of the central 
node (filled circle) share another vertex (dotted empty circle) 
as nearest neighbor. 



Fig. shows what we shah cah a tertiary quadrilat- 
eral: the three nearest neighbors of the central node are 
all connected to another common node. Obviously, the 
presence of a tertiary quadrilateral implies Q2 — (2) ~ ^ 
secondary quadrilaterals. In a fourth-order quadrilateral 
(Fig. 13)), 4 nodes share two common vertices as nearest 
neighbours, implying the presence of Q3 — {^^ — A ter- 
tiary quadrilaterals and Q2 = (2) = ^ secondary quadri- 
laterals. In a network of degree K > 4, quadrilaterals up 
to order K can in principle be found. 

A one-dimensional lattice with additional connections 
between the nodes i and i -t- 3 for all i (instead of i -f 2 as 
in the ring investigated up to now, Fig.QJi) possesses the 
neighbourhood structure shown in Fig. llUI i.e. it is char- 
acterized hy E — Qi — 0,Q2 — 8,Q3 — 2 and Q4 — 0. 
By applying our formalism to this case and to a network 
that has fourth-order quadrilaterals, Eq. @ generalizes 
to 



1 = A 



16A- (6 + 2S + Q1 + Q2)A2 

+ (4 + Qi + Q3)A3-(1 + Q4)A^ 



(10) 



the coefficients of order 3 and 4 being modified only. 



C. Introducing disorder 

The networks considered up to now lack in the small 
world phenomenon, a property characterizing social net- 
works on which the epidemic process is occuring. By 
starting with a ring-like network (Fig. ^) and repeating 
the rewiring procedure described in section II a certain 
number of times, we obtain graphs of fixed degree K ^ 4 
that are simultaneously highly clustered, and in which 
the average distance between pairs of nodes is small 3] . 

The left part of Fig. ^] reports the simulation result 
for the equilibrium prevalence of the epidemic process on 
the disordered ring. Systems of size N — 10'* were used, 
and the rewiring procedure was repeated n = 100 times. 
For completeness, the two limiting cases (fully ordered 
ring and random network) are also depicted. This panel 
shows the considerable effect on the steady-state spread- 
ing behavior of the rather small number of rewirings. 




FIG. 10: This figure visualizes how the nearest (empty solid 
circles) and second neighbors (empty dotted circles, vertices 
{-3,-1,1,3}) of an arbitrarily chosen node (filled circle, node 0) 
are arranged in a one- dimensional lattice with additional con- 
nections between sites 3 units apart. The sets {0, {-1,1, 3}, 2} 
as well as {0,{-3,-l,l},-2} are forming tertiary quadrilaterals. 



The right panel of this figure depicts the one-step site 
approximation (predicting Ac = 0.25 for all cases) and 
the numerical solutions of the two-step description ||SJ. 
It has to be noted that for the partially rewired ring lack- 
ing in strict homogeneity, Eq. (0) was solved at every 
node, therefore involving the set Pi{xi), i — 1,2, ...,7V, 
the resulting prevalence being given by ^/Nj^iLi ^i(l)- 
Again at the double-step level, the networks in question 
are treated differently, as it could already be observed 
in Fig. 121 Here the curve corresponding to the partially 
randomized ring lies closer to the original network in pro- 
portion to the simulation result. This is due to the small 
world property which can have a considerable effect on 
the location of the onset of the epidemic. Obviously, the 
simulation result uncovers the real effect of this global 
topological property whereas at the two-step level, it is 
the slightly poorer loop structure that accounts for the 
corresponding shift in the epidemic threshold. 



Of course the quantities -E, Qi, (32, Q3 and Q/^ are no- 
longer reasonable for a partially randomized network due 
to the lack of strict homogeneity, but rather its local or- 
dering properties can be quantified by averaging these 
values over the entire network. The emerging topolog- 
ical parameters E and Qi {i = 1,2,3,4) are essentially 
the clustering- (up to the factor K{K — l)/2 — 6) and 
grid- coefficients, i.e. the densities of triangles and loops 
of length 4 Q, ■ We may therefore replace E and the 
number of quadrilaterals (of the different orders) in Eq. 
(|10|l by its mean-values, yielding the following estimate 



0.2 




0.35 



FIG. 11: Steady-state prevalence for the ring-type network having undergone different degrees of randomization. Left: Simu- 
lation result for the original network (o), its partially rewired version (A) and the entirely random network (★). Right: The 
single-step site approximation (thin solid) treats all the networks in question identically whereas at the two-step level, the 
MoUoy Reed network (dotted), the partially randomized ring (solid) and the fully ordered ring (dashed) appear in the same 
sequence as in the left panel. 



for the epidemic threshold condition 



1 = A 



+ (4 + Qi + Q3)A^-(l + Q4)A^ 



(11) 



For our partially randomized ring, we have E — 
2.883, Qi = 1.886, Q2 = 1-958 and Q3 = Q4 = 0, leading 
to Ac — 0.2892. This value corresponds approximately to 
where the corresponding curve in the right part of Fig. 
II II meets the x-axis. 



V. ARBITRARY DEGREE 

The implications of our two-step description have been 
illustrated for homogeneous networks of degree 4 in the 
previous section. This was a convenient choice as there 
exists a number of familiar simple graphs obeying P{k) — 
Ski, differently ordered. Of course our formalism enables 
us to generalize the obtained threshold condition I|1U|) 
to an arbitrary degree K, which is the subject of this 
section. 

Let us again look at a fully treelike network, using 
Eqs. ijnji and © as guidelines. The A^-coefficient 16 in- 



corporating the degree distribution is simply 4 • 4 since 
(/„) = 4P + 0(p2) and a runs from 1 to 4. The re- 
maining coefficients -6, 4 and -1 correspond to the bino- 
mial coefficients — (2) , (3) and — (^) . Indeed the thresh- 
old equation for a treelike network of degree K derived 
by Eq. ® is 



1 = A 



K 



Aif 2 - ^ A« 



K 



(12) 



Repeating the graph developments for homogeneous 
networks characterized by different values of K and vary- 
ing loop structures reveals that the very same correction 
terms enter into Eq. ifT^ . yielding 



1 = A|e(if 

-e{K 

+e(if-3) 

K 



. .K 

K 
3 

K 



2E + Q1+ Q2 

Oi + OajA^ 



(13) 



i-xr} 
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where E is again the number of connections between the 
nearest neighbors of an arbitrarily chosen node, Q„ de- 
notes the number of quadrilaterals of order n and Q{x) 
is the step-function defined by 




1 if a; > 
otherwise. 



The equation that describes the role of local ordering in 
a disordered network subject to P{k) = 5k. k is again ob- 
tained simply by replacing the corresponding quantities 
by its mean-values [E ^ E, Qi ^ Qi) in Eq. ((T^ . pro- 
viding an improved estimate for the epidemic threshold. 

VI. CONCLUSION 

The spreading of an infectious disease was modeled as a 
dynamical process on top of a contact network. We used 
a discrete-time version of the simple SIS-model, that is 
infected nodes recover with probability At, and suscep- 
tible nodes become infected with probability AAi if they 
are connected to at least one infected nearest neighbor. 
As far as the connectivity patterns underlying the popu- 
lation are concerned, we chose homogeneous networks as 
starting point and introduced an arbitrary degree of dis- 
order by an appropriate rewiring procedure not affecting 
the degree distribution P{k) = Sk,K- 

Describing the epidemic dynamics of the entire popu- 
lation as a Markovian process, we derived a two-step de- 
scription that takes into account temporal correlations. 
This approach revealed to be very prolific if one wants 
to unravel the role of loops of short length in the contact 
network regarding epidemic spreading. Indeed it leads 
to a subgraph development where the complete graph 
involves the connectivity patterns of two hierarchies of 
nearest neighbors (of an arbitrarily chosen node). Within 



this novel approach serving to tract a probabilistic sys- 
tem, the local topology, be it treelike or be loops of length 
3 or 4 present, therefore enters very naturally. The ana- 
lytically obtained condition for the location of the onset 
of the epidemic then serves as a guiding equation elu- 
cidating the role of clustering and grid-like ordering in 
epidemic spreading. 

In principle it is possible to apply our two-step de- 
scription to more complex networks where different de- 
grees are present, uncovering the effect of the degree- 
dependent densities of triangles and loops quadrilaterals 
on the critical value. Likewise, loops of length up to 2n 
are expected to enter within an n-timestep description, 
also providing a natural classification of them. However 
the major insight gained by the strategy of exploring tem- 
poral correlations is best illustrated as it was done in this 
paper. 
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APPENDIX: FULL SUBGRAPH 
DEVELOPMENTS 

Here we show the complete subgraph developments for 
the square lattice (Tab. IIII|I . the Kagome lattice (Tab. 
IIV(I and the ring-type network of Fig.ni fTab. lV)l . Every 
subgraph corresponds to a term in Eq. (jSJ, its contribu- 
tion is obtained by the procedure illustrated in section 
IVA. The A • A"-coeffLcient of the threshold equation is 
obtained by summing all the 0{P) contributions (tak- 
ing into account the multiplicities) of the n-th-order sub- 
graphs. 
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TABLE III: Subgraph development for the square lattice. All the terms in Eq. Q are symbolized by a specific subgraph, its 
order being given by the number of filled circles. The A • A^-coefficient -10, as an example, is obtained by summing the various 
0{P) contributions, that is2-4 + l- 2 + 0- 4 = 10, and the negative sign comes from Eq. @. The A ■ A^- and A ■ A'*-coefEcients 
are unaltered with respect to the treelike topology although other subgraphs enter into the development. 



order (n) 



subgraph 



contribution 
multiplicity 



4P + 0{P'^) 
4 



2P + 0(p2 
4 



P + 0{P^ 

4 



P + 0{P^ 
1 



subgraph 



contribution 
multiplicity 



P + 0{P^) 
2 



2P2 + 0{P^) 
4 



(B (I G) 



e ® o 



^P^ + 0{P^) 
4 



subgraph 



contribution 
multiplicity 



Q fl) 



(3- ffi 



9P^ + C'(P^) 
4 



o ffl ffl 



ffl (t) o 



18P^ + 0{P^) 
2 



A ■ A"-coeff. 



16 



-10 
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TABLE IV: The subgraphs of aU the orders for the Kagome lattice. See Tab. Illll for how the coefficients are obtained and as 
far as further details are concerned. 



order (n) 


I 


2 


3 


4 


subgraph 

contribution 
multiplicity 






G V 




G \- — G 








— i*- — © 


G Hl(- -llh G 


4P + C'(p2) 
4 


P + 0{P^) 
6 




P + 0{P^) 

4 


P + 0{P^) 
1 


subgraph 

contribution 
multiplicity 




Q 






G -^1^^ 0------ 






^^^^^ 


Q -fh €!; 

P + 0{P^) 
4 


5P^ + 0{P^) 
4 


SP"* + C(P^) 
4 


subgraph 

contribution 
multiplicity 












G ^» 

8P^ + ( 
2 


— r — ® 

«; 

9(P'*) 


A ■ A"-cocff. 


16 


-10 


4 


-1 
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TABLE V: The full subgraph development for the ring- type network. See Tab. IIIl| for the derivation of the A • A"-coefiicients. 
With respect to the two lattices treated above, the A ■ A^- and A • A'^-coefficients are -16 and 6. 



order (n) 



subgraph 



contribution 
multiplicity 







subgraph 



contribution 
multiplicity 






subgraph 



contribution 
multiplicity 






subgraph 



contribution 
multiplicity 




3 




A ■ A"-coeff. 



16 



-16 



